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Abstract: We implement the Zee-Babu model for the neutrino masses and mixings by 
incorporating a scalar dark matter X. The singly and doubly charged scalars that are new 
in the Zee-Babu model can explain the large annihilation cross section of a dark matter 
pair into two photons as hinted by the recent analysis of the Fermi 7-ray space telescope 
data, if new charged scalars are relatively light and have large couplings to a pair of dark 
matter particles. These new scalars can also enhance the B(H — > 77), as the recent LHC 
results may suggest. The dark matter relic density can be explained. The direct detection 
rate of the dark matter is predicted to be about one order of magnitude down from the 
current experimental bound. However, it turns out that neutrino masses are too small 
within the parameter space fitting the 130 GeV 7-ray excess. There should be additional 
contributions to the neutrino masses and mixings without new extra charged particles, and 
the Type-I seesaw can do this job. 
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1 Introduction 

Although it is well known that the dark matter (DM) constitutes about 23% of the total 
mass density of the universe, i.e. f^DM^ 2 = 0.1123 ±0.0035 [1], its existence has only been 
inferred from the gravitational interaction. And its nature is still unknown. If the DM 
is weakly interacting massive particle (WIMP), it may reveal itself via non-gravitational 
interactions, for example, by pair annihilation into ordinary standard model (SM) particles 
including photon [2]. In this case, the relic DM abundance is roughly related to the pair 
annihilation cross section at freezeout, (cu)th) as 

Recently Refs. [3, 4] claim that the Fermi 7-ray space telescope may have seen excess of 
the photons with ~ 130 GeV from the center of the Milky Way compared with the 
background. Interpreting its origin as the annihilation of a pair of DM particles, they 
could obtain the annihilation cross section to be about 4% of that at freezeout: 

(<™) 77 = 0.042(o-v) th = 0.042 pbc. (1.2) 

Since DM is electrically neutral, the pair annihilation process into photons occurs through 
loop-induced diagrams. Naively we expect 

M 77 _ ^« em x2^ io _ 5 _ 
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So the observed value in (1.2) is rather large, and we may need new electrically charged par- 
ticles running inside the loop beyond the SM. Many new physics scenarios were speculated 
within various CDM models by this observation [5]. 

The so-called 'Zee-Babu model' [7-9] provides new charged scalars, h + ,k ++ at elec- 
troweak scale, in addition to the SM particles. These new charged scalars carry two units 
of lepton number and can generate Majorana neutrino masses via two-loop diagrams. The 
diagrams are finite and calculable. The neutrino masses are naturally small without the 
need to introduce the right-handed neutrinos for seesaw mechanism. One of the neutrinos 
is predicted to be massless in this model. Both normal and inverse hierarchical pattern 
of neutrino masses are allowed. The observed mixing pattern can also be accommodated. 
The model parameters are strongly constrained by the neutrino mass and mixing data, the 
radiative muon decay, \x — > cy, and r — > 3[jl decay [10]. 

It would be very interesting to see if the new charged particles in the Zee-Babu model 
can participate in some other processes in a sector independent of neutrinos. In this 
paper, we minimally extend the Zee-Babu model to incorporate the DM (See also [11]). 
We introduce a scalar dark matter X with a discrete Z2 symmetry under which the dark 
matter transforms as X — > —X in order to guarantee its stability. The renormalizable 
interactions between the scalar DM X with the Higgs field and the Zee-Babu scalar fields 
provide a Higgs portal between the SM sector and the DM. We show that the Zee-Babu 
scalars and their interactions with the DM particle can explain the gamma-ray excess as 
well as the DM relic density. The branching ratio of Higgs to two photons, B(H — > 77), 
can also be enhanced as implied by the recent LHC results [12]. The spin-independent 
cross section of the dark matter scattering off the proton, a p , is less than about 1 x 10 -9 
pb, which can be probed at next generation searches. 

This paper is organized as follows. In Section 2, we define our model by including 
the scalar DM in the Zee-Babu model, and consider theoretical constraints on the scalar 
potential. In Section 3, we study various DM phenomenology. We calculate the dark matter 
relic density and the annihilation cross section (crf) 77 in our model. We also predict the 
cross section for the DM and proton scattering. We also predict the branching ratio for 
the Higgs decay into two photons, B(H — > 77). And we consider the implication on the 
neutrino sector. We conclude in Section 4. 



2 The model and theoretical constraints on the scalar potential 

2.1 Implementation of Zee-Babu model with scalar dark matter 

We implement the Zee-Babu model for radiative generation of neutrino masses and mixings, 
by including the scalar DM X with Z2 symmetry X — > —X. All the possible renormalizable 
interactions involving the scalar fields are given by 

£ = ^Babu + -Cffiggs+DM (2.1) 
£ B abu = fabl T atCl{ L e l3 h + + h' ab ll R Cl bR k ++ + H.C. (2.2) 

-£ffi gg s + DM = Xh(H^H) 2 + \\xX A + X h (h + h-) 2 + \ k {k ++ k-f 



-2- 



+ ^X H xH^HX 2 + X Hh H^Hh + hT + X Hk H^Hk ++ k- 

+ ±X Xh X 2 h+h~ + h Xk X 2 k ++ k- + X hk h+h'k ++ k- 
+ {pihkh + h + k~ + h.c) 

- ii 2 H rfH + ^ 2 X X 2 + fih+hr + ^k ++ k-. (2.3) 

Note that our model is similar to the recent model proposed by J. Cline [6]. However 
we included the interaction between the new charged scalar and the SM leptons that are 
allowed by gauge symmetry, and thus the new charged scalar bosons are not stable and 
cause no problem. 

The original Zee-Babu model was focused on the neutrino physics, and the operators of 
Higgs portal types were not discussed properly. It is clear that those Higgs portal operators 
we include in the 2nd line of (2.3) can enhance H — > 77, without touching any other decay 
rates of the SM Higgs boson, as long as and k^ are heavy enough that the SM Higgs 
decays into these new scalar bosons are kinematically forbidden. 

2.2 Constraints on the potential 

We require n 2 x , n\ and //| to be positive. Otherwise the imposed Z2 symmetry X — > —X 
or the electromagnetic U(l) symmetry could be spontaneously broken down. Since the 
masses of X, k ++ and h + have contributions from the electroweak symmetry breaking as 

m x = (4c + ^W, 

m l+ = Vh + \ x HhV 2 , 

= A + \^HkV 2 , (2.4) 

we obtain the conditions on the quartic couplings 

2m 2 x ^ m \+ ^ m 1++ 
Xhx < — 5—, X H h < — n— , X Hk < 5 — • (2-5) 

yZ yA yZ 

We note that the above conditions are automatically satisfied if the couplings takes negative 
values. In such a case, however, we also need to worry about the behavior of the Higgs 
potential for large field values. For example, if we consider only the neutral Higgs field, 
H , and the dark matter field, X, we get 

V ~ ^X H H 4 + ^X X X 4 + ^X HX H 2 X 2 , 




Xh \Xhx 
\Xhx Xx 




(2.6) 



for large field values of H and X. If the potential is to be bounded from below, every 
eigenvalue of the square matrix of the couplings in (2.6) should be positive, whose condition 



x We use the same notation with the Higgs doublet. 
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is 



(2.7) 



This means that even if Xhx is negative, its absolute value should not be arbitrarily large 
because Xh = m? H /2v 2 fa 0.13 (m# « 125 GeV) and Xx is bounded from above so as not 
to generate the Landau pole. For example, the renormalization group running equations 
(RGEs) of Xh, Xx and Xhx are given by 



where the dots represents other contributions which are not important in the discussion. 
The complete forms of the /3-functions of the quartic couplings are listed in Appendix A. 

The approximate solution for Ax in (2.8) shows that the Landau pole is generated at 
the scale Q = Qew exp (l//3ff Ax(Qew)) {Ph = 18/167T 2 ). If we take the electroweak scale 
value of the Higgs quartic coupling to be Xx(Qew) ~ 5, the cut-off scale should be around 
1 TeV. The general condition for the bounded-from-below potential for large field values 
is that all the eigenvalues of the matrix 



should be positive. 

In the following discussion, we require that all scalar quartic couplings (Aj) be per- 
turbative up to some scale Q. To this end, we solve the one- loop RGEs of those quartic 
couplings given in Appendix A. For the moment, we do not include new Yukawa couplings 
defined in Eq. (2.2), and we adopt the criterion Aj(Q) < 47T in this analysis. In Fig. 1, the 
perturbativity bounds are shown in the X X h(k)-^Hh(k) plane. We take Q = 1, 3, 10 and 15 
TeV, which are denoted by the red curves from top to bottom. For other parameters, we 
fix X Hh = X H k, ^Xh = Axjfe, X hk = X H x = and X x = X H (~ 0.13). As explained above, a 
certain negative value of Xn k {h) mav cause the instability of the Higgs potential. To avoid 



this, we set A^ = X 2 Hh /(2X H ) and X k = X Hk /(2X H ) for X H k{h) < 0. For X Hk ^ > 0, on the 
other hand, Xh = X k = Xh is taken. As we see from the plot, Xxk(h) — 7 — 11 is possible if 
Q = 1 TeV. 

The theoretical arguments (2.5) and (2.7) restrict Xhx to lie roughly to the range, 




(2.8) 



( Xh \XhX ^Hh ^Hk ^ 

\^hx Ax Xxh Axfe 
Atf/i Xxh 4Ah 2Xhk 
\ ^Hk Ax/i 2Xhk 4Afc J 



(2.9) 



(—1.6,0.6). Similarly, we have Xnh(k) ^ 0.7 for m h +n i ++\ = 150 GeV. 
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Figure 1. The perturbativity bounds, Xi(Q) < Att are shown. The each curve denotes Q = 1, 
3, 10 and 15 TeV from top to bottom. We take Xhh = \Hk, Xxh = Xxk, X^k = Xhx = and 
Xx = Xh{— 0.13). For the negative Xhu, we set Xh = X 2 Hh / (2Xr) and Xk = X 2 Hk /{2Xij) while 
Xh = A fc = X H for positive X H k- 



3 Dark matter phenomenology and H — y 77 

3.1 XX 77 and Fermi/LAT 130 GeV 7-ray excess 

The annihilation cross section for XX — > 77 is given by 

_£|M| 2 



(av) 



77 



6Airm x ' 



where the amplitude-squared summed over the photon polarization is 

„2 



+ 



a 

Am x — m 2 H 



^XhTh+A (T h +) - AXxkT k ++A Q (T k ++) 



Q 2 / 

y —T W (Q 2 N c A 1/2 {T t ) + Ai(tw) 



with Tj = m? x /m 2 {i = h + , k ++ ,t, W). The loop functions are 



^1/2 (t) 



A (r) = - r-/(r) 



r+(r-l)/(r) 



2r 2 + 3r + 3(2r- 1)/(t) 



(3.1) 



(3.2) 



(3.3) 
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(a) A m =^uh -^Hk -Q (b) Ahx^Ahh =^Hk =0-33 




100 110 120 130 140 150 100 110 120 130 140 150 



m /j+ (GeV) m h+ (GeV) 

Figure 2. Contour plot of (<rv) 17 = (2, 1, 0.5, 0.2) x 10~ 27 cm 3 /s (from above) in (m^+, Xxh) plane. 
We set ra x = 130 GeV, m H = 125 GeV, m k = 500 GeV and X xk = 5, Xhx = X Hh = X Hk = 0(0.33) 
in the left (right) panel. 



where 



f(r) 



arcsin 2 \pr. (t < 1] 



l0 § ^^Epr -^ , (r>l). 



(3.4) 



Although the contribution of the doubly-charged Higgs k ++ to (<jv}y~ ( is 2 4 = 16 times 
larger than that of the singly-charged Higgs h + when their masses are similar to each 
other, this option is ruled out by the recent LHC searches for the doubly-charged Higgs 
boson [13]. Depending on the decay channels, the 95% CL lower limit on the mass of 
the doubly-charged Higgs boson is in the range, 204-459 GeV. To be conservative, we set 
m/j++ = 500 GeV. In Figure 2, we show a contour plot for the annihilation cross section 
into two photons: (<tu) 77 ~ (2,1,0.5,0.2) x 10 _27 cm 3 /s (from above) in the (m h +, Xxh) 
plane. We set rax = 130 GeV, ran = 125 GeV, m k ++ = 500 GeV and Xxk = 5, Xhx = 
Xhh = Xuk = (0.33) in the left (right) panel. We can see that even with Xhx = 0.33, the 
process XX — )■ H — > 77 can be significant. However, this value of Xhx produces too large 
cross section to explain the current CDM relic density, and and we would need another 
DM components. We can see that the Fermi result can be explained when lies in the 
range 100-135 (150) GeV for Xxk ~ 5 with the cross section (cru) 77 w 1(0.2) x 10~ 27 cm 3 /s. 
In this case the cut-off scale can be as large as 15 TeV (See Fig. 1). This mass range is 
within the reach of the on-going LHC experiments. And the model can be tested soon. 

3.2 Direct detection rate and thermal relic density 

Contrary to J. Cline's model [6], the DM relic density in our model is not necessarily 
correlated with the {av ) 77 , since it is mainly determined by Xhx for relatively heavy scalars 
(> 150 GeV). In this case the main DM annihilation channels are XX — > H — > SM particles, 



-6- 



m h =l25 GeV,m x =130 GeV 




-10 -5 5 10 

Figure 3. The contour plot of S] D m'* 2 = 0.1123 (red lines) and {av) 11 = 0.2 x 10~ 27 cm 3 /s 
(black lines) in the (Xxii-.Xhx) plane for the choices m h + 150, 140, 130 GeV (solid, dashed, dotted 
lines). For other parameters we set mx = 130 GeV, m# = 125 GeV, mk = 500 GeV, Xxk = 5, 
Xflh = Xffk = 0.5. 

where the SM particles are W + W~, ZZ, bb, etc. As m h +t k ++\ becomes comparable with 
mx, the XX — > h + h~(k ++ k ) modes can open, even in cases mx < m h +(m k ++) thanks 
to the kinetic energy of X at freeze-out time. This can be seen in Figure 3, where we show 
the contour plot of Odm^ 2 = 0.1123 (red lines) in the (Xxh^hx) plane for the choices 
m h + = 150, 140, 130 GeV (shown in solid, dashed, dotted lines respectively). We fixed 
other parameters to be mx = 130 GeV, m# = 125 GeV, m k = 500 GeV, Xxk = 5, 
^Hh = ^Hk = 0.5. For m h + = 130 GeV, the annihilation mode XX — > h + hr dominates 
even for very small coupling Xxk (the red dotted line). The black vertical lines are the 
constant contour lines of (crv)yy = 0.2 x 10 _27 cm 3 /s. Again the solid (dashed, dotted) line 
is for m h+ = 150 (140, 130) GeV. 

For m h + = 150 GeV, black line meets the red line, which means that if (<rf) 77 turns 
out to be about 0.2 x 10 _27 cm 3 /s in future experiments, our model can both explain the 
Fermi/LAT result and saturate the relic density. Otherwise, the (<7t>) 77 remain as large as 
the current value, the relic density is too small to explain the observed value and we need 
other DM components. For m^+ = 130,140 GeV, even (o"^) 77 = 0.2 x 10 _27 cm 3 /s is too 
large to saturate the relic density. 

Figure 4 shows the cross section of dark matter scattering off proton, a p , as a function 
of Xhx (red solid line) and a p = 1.5 x 10 -8 pb line (black dashed line) above which is 
excluded by XENON100 [14] at 90% C.L. This cross section is determined basically only 
by Xhx at tree level by the SM Higgs exchange, when we fix mx = 130 GeV. We can see 
that Xhx ^5 0.1 to satisfy the XENON100 upper bound. We see that the cross section is 
close to the current bound and can be probed by near- future experiments when Xhx ^ 0.02. 
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Figure 4. The spin- independent cross section of dark matter scattering off proton, a p , as a function 
of Xhx (red solid line) and a p = 1 X 10 -8 pb line above which is excluded by XENON100 (black 
dashed line). We take mjj = 125 GeV and mx = 130 GeV. 



If the new scalars are light and dominate the DM annihilation cross section for the relic 
density, we have Xhx ^ 0.02 and the a p is one order of magnitude less than the current 
bound. 



3.3 H 



77 



In our model the decay width of H 
are intact: 



77 is modified, whereas other Higgs decay widths 



F(H -> 77) 



GFa 2 m H 



+ 



128\/27r 3 
^Hhmwv 



]T Q}N c A 1/2 (r f ) + A^tw) 

f=t,b 



9 m l 



A\ Hk m w v 

2 A 0{Th+) + —Z2 A 0{T k ++) 



9 m l+- 



(3.5) 



where r» =rn H / Amj (i = f,W,h + ,k ++ ). 

In Figure 5, we show contour plots for constant T(H — > r yy)/T(H — > 77) SM (black 
solid lines) and T(H — > Zj)/T(H — > Zj) SM (black dashed lines) in the (\fih,^Hk) plane. 
For this plot we set m/j+ = 130 (150) GeV for the left (right) panel and fixed m k ++ = 500 
GeV. The shaded regions are disfavored by (2.5) (blue) and by (2.7) (yellow). Given that 
^hx ^ 0.1 by the direct detection experiments, the ratio depends basically only on the 
coupling constants and Xhu as well as the masses m h + and m k ++. In our model the 
ratio is not necessarily correlated with the (crv) 77 which are controlled by Xxh and Xxk- 
We can conclude that 

0.54 < T(H -> Tf)/T(H -> 7 7) SM < 1-45 (1.35) 
0.91 < T(H Zj)/T(H -»• Z 7 ) SM < 1.11 (1.08) 
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Figure 5. A contour plot for constant T(H — > 77)/r(iJ — > 77)™ (black solid lines) and T(H — > 
Z"f)/T{H —> Z^f) SM (black dashed lines) in the (\Hh,^Hk) plane. The shaded regions are disfavored 
by (2.5) (blue) and by (2.7) (yellow). We set m h + = 130 (150) GeV for the left (right) panel and 
fixed m k ++ = 500 GeV. 



for the left (right) panel. That is, the H — > 77 channel can be enhanced (reduced) signifi- 
cantly, whereas the H — > Z7 channel can change only upto ~ 10%. 

3.4 Implications for neutrino physics 

So far, we did not consider the constraints from neutrino sector and charged lepton flavor 
violation. In fact, these constraints are rather severe, if we assume that the observed 
neutrino masses and mixings are entirely from the Zee-Babu mechanism. One cannot 
afford light or k^, because of the constraints from charged LFV: m h + > 240 GeV 
from B([i — > e^f) < 2.4 x 10 -12 [15] and m^++ > 770 GeV from the upper bound on r — > 3/U 
decay [10]. 

If we insist that the Fermi/LAT excess is due to the light h^- loop, then we get fihk ^ 14 
TeV, which is inconsistent with the constraint /i/^ < 450 GeV from the vacuum stability 
bound [9]. And there should be additional contributions to neutrino masses, such as from 
dim-5 Weinberg operators. If these dim-5 operators are induced through Type-I seesaw 
mechanism, the new physics would not affect our conclusion. On the other hand, if the 
dim-5 operators are induced by TeV scale seesaw, then the new physics from TeV scale 
seesaw might affect our conclusions. 

Although it is not very satisfactory that the original Zee-Babu model with scalar dark 
matter cannot explain both the Fermi/LAT 130 GeV 7 ray excess and neutrino physics 
simultaneously, it would be more natural to consider the Zee-Babu model as a low energy 
effective theory. Then it would be natural there could be new contributions to the neutrino 
masses and mixings from dim-5 operators. The only relevant question would be whether 
those new physics would affect the Fermi/LAT 7-ray or not. If the new physics is Type-I 
seesaw, there would be no new charged particles so that our conclusion would remain valid. 
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4 Conclusions 



We have minimally extended the 'Zee-Babu model' [7-9] incorporating a scalar dark matter 
X with Zi symmetry under X — > — X. If the scalar dark matter X has a mass around 
130 GeV, it can explain a recent claim that Fermi gamma-ray space telescope has seen an 
excess in gamma line at 130 GeV. If the excess is interpreted as the annihilation of a dark 
matter pair into two photons, the estimated cross section (av)~~ ~ 1 x 10~ 27 cm 3 /s can be 
explained by the contribution of the singly- and/or doubly-charged Higgs present in the 
model. 

We have shown that the present constraint on the couplings Xxk and Xxh which mix 
the dark matter and charged Higgs is not so strong and they can enhance the annihilation 
cross section of XX — > 77 large enough to accommodate the recent hint. On the other hand 
the couplings which involve the SM Higgs H are strongly constrained by the theoretical 
considerations in the Higgs potential and the observations of dark matter relic density and 
dark matter direct detections. The strongest upper bound on the Xhx coupling is about 
0.1 which comes from the dark matter direct detection experiments. For the \Hh^Hk 
which mix the SM Higgs and the new charged Higgs, the theoretical bound becomes more 
important. If we require the absolute stability of the dark matter by the Z2 symmetry 
X — > — X and the absence of charge breaking, we get the upper bound of \Hh, ^Hk to be 
about 0.7 for the charged Higgs mass around 150 GeV. To evade the unbounded-from-below 
Higgs potential we need to have Xjih, ^Hk ^ —1.6. 

With these constraints the B(H — > 7(^)7) can be enhanced up to 1.5 (1.1) or sup- 
pressed down to 0.5 (0.9) with respect to that in the SM. The direct detection cross section 
of X is suppressed by more than an order of magnitude to the current XENON100 bound 
if we want to saturate the relic density. 

However the neutrino sector cannot be described by the Zee-Babu model only, and there 
should be additional contributions to the neutrino masses and mixings such as dimension-5 
Weinberg operator from type-I seesaw mechanism. As long as the new mechanism does not 
involve new charged particles around EW scale, our predictions for 130 GeV 7-ray excess 
and the enhanced B(H — > 77) would be still valid. 



Acknowledgment s 

This work is partly supported by NRF Research Grant 2012R1A2A1A01006053 (PK, SB). 



A One-loop /3 functions of the quartic couplings 

Here, we give the renormalization group equation and the one-loop /3 functions of the 
quartic couplings: 
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with 
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^Xjih^Hk + 8Xhk(hi + Afc) + Xxh^Xk + ^8gf — 30Xhkgl 



(A.2) 
(A.3) 

(A.4) 
(A.5) 

(A.6) 

(A.7) 
(A.8) 



12X H X HX + 2X Hh X X h + 2A fffc Axfc + 6AxAhx - ^hx{ ^(3#1 + # 2 ) - 6y 2 
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